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' Abstract 

' Network (as a general notion) is not a mathematical object - there is no even any definition. However, 

there is a lot of good rigorous mathematics for well-defined classes of networks. In sections 1-3 we give a 
short overview of classes of networks which interested the authors for some time. In section 4 we consider 
' ' in detail a new class of networks, related to markets with many agents. 

1 Random field dynamics on a fixed graph 

■ The basic element of most networks is a graph G with the set V — V{G) of vertices and the set L = L{G) of 
\ links (lines, edges). Second basic element is a function s — f{v) : V ^ S with values in some space S The 

^ ■ elements of S may be called marks, spins, field values, queues etc. The function / is subjected to random 
dynamics. 

I ■ Simplest example is (an earlier stuff) random walks on graph, where f — everywhere except one point 

^ I where the particle is situated. This is related to electric networks, see for example [6J. 

00 ' In general there are two different situations. First one is a local continuous time Markov dynamics given by 

On I infinitesimal transitions. Classical reference is [9], mostly such processes model stochastic dynamics of particles 

■ or spins. The latter are related to Gibbs random fields (invariant measures for this dynamics) on graphs, see 
' Hn] and references therein. 

(n: 

, Queuing, communication and transportation networks The simplest case is when the particles jump 

' (from one node to another) freely without seeing each other, the only interaction is only through queues at the 

IL" , nodes, where they spend some time. There are two main theories concerning such class of networks: 
^ , 

• 1. Most popular - Jackson network (1963) and its generalizations (Gordon-Newel, BCMP). This theory gives 

I explicit formulas for the stationary distribution and is the origin of many other analytical results. One of 

' the applications is to describe jams and phase transitions in communication |11) and transportation |12) 

networks. 

2. Stability theory (1968-1995) exhibits in many cases of large time qualitative behavior. If the walking 
clients are identical then it is described by random walks in orthants and strongly uses It uses Lyapounov 
functions, Euler scaling (fluid approximation), ergodic theory of dynamical systems and Lyapounov ex- 
ponents, see [13|. If the walking clients can be of finite number of types then the corresponding theory [4] 
is the union of the one type case and the theory of random grammars (see below). 

For more sophisticated restrictions - network protocols (TCP etc.) - there are many partial results but no 
comparable (deeply elaborated) mathematical theories. 
Chemical kinetics - mean field Markov chain 

Mean field network means that there is no specified local structure on the graph. Example of such theory 
is the chemical kinetics. It describes the following situation. Molecular types are indexed by V , riy - number of 
molecules of type v 

fii + ... + n\v\ = N 

There are also reaction types r = 1,2, R, formally - multigraph defined by finite number of equations 

S^rMy = 
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where My - molecule of type w, s^r - stoihiometric coefficients of molecule type v in reaction of type r , negative 
for substrates, positive for products. Reaction rates (continuous time Markov chain) are given by 
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D:S„r<0 



for the jump (transition) 
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To get ODE of classical chemical kinetics 



lit 



'V 



Qvr (ci ; ■ • ■ 1 



c\v\) 



r 



for some polynomials Qr, in the limit N ^ oo 
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one uses canonical scaling of reaction rates 
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v:s^r<0 



To deduce chemical thermodynamics is more difficult [14], one should, together with molecular types v, 
introduce more degrees of freedom: kinetic energy Ty^i and internal energy Ky i of j-th molecule of type v. Also, 
one should define more complicated mean field dynamics - introduce energy mechanism in reactions. As there is 
kinetic energy - there should be Newtonian movement, and the dynamics become mixed: local + mean field. 
Molecule move freely (as in ideal gas) but kinetic energies randomly interchange with internal energies. 

Network homeostasis [15] Network is defined by 

1. large graph G of compartments, this graph G has metrics and the boundary, 

2. in any compartment chemical kinetics is defined, that is there are molecules with chemical reactions, 

3. there is transport of molecules between compartments 

4. there is input and output of molecules on the boundary 

Under some conditions (the main is that reactions are unary) it is possible to prove that far from the boundary 
there is equilibrium - concentrations almost do not change with change of input. 

2 Dynamics of graphs and of marked graph 

Earlier the science of random graphs considered mainly the properties of graphs with fixed number of vertices 
and/or random number (for example Bernoulli) of edges, see for example [6l [TBI [191 [8] . The simplest dynamics 
(appending edge by edge) appeared already in 0, see also [3]. What more general dynamics on graphs one should 
study? First of all, it is more reasonable to consider evolution of marked graphs. Most general dynamics of 
marked graphs (local random dynamics of a graph, jointly with a field on it) is called random graph grammars 
[21 [31 [T] and [3T1 [52] . It appears to be quite natural in connection with the emerging new physical theories 
|16[ [TTI [25] and social networks [20] . Namely, if eventually the local space-time appears to be discrete, then the 
most natural language for it is a graph with some physical fields on it. The dynamics of the space time is local. 
The example is the following. 
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Macrodimension of a graph - invariant of local dynamics We consider infinite (countable) graphs G. 
Let On{v) be the neighborhood of vertex v of radius n. Put 

Dr,{v) = i^ii^M^, = lim sup Dn{v),D{v) = lim inf Dn(v) 

inn n— i-oo n— s-oo 

If for all V 

D{v) = D{v) = Ds 

then Ds is called scaling macrodimension of graph G. For example any homogeneous lattice in euclidean 
space R'^ has scaling macrodimension Ds — d. Note that there are many other definitions of variants of 
macrodimension: connectivity, Hausdorf, entropy, inductive macrodimension. 

Denote Gm the class of connected graphs where each vertex has degree < M . Let U - any local dynamics 
(graph grammar). 

There is the following result [IB] . If for some sufficiently large M U leaves the class Gj^j invariant and the 
corresponding Markov chain is locally reversible then the scaling macrodimension is an invariant. 
Local reversibility means that Kolmogorov cycle criteria relations 

— 1 — hll. 

follow from such relations of bounded length. 

Random graph grammars Consider words a — xi...xn (ordered sequences of symbols), where xn belongs 
to some finite alphabet A. Grammar is defined by the list Sub of productions (allowed substitution types) 

Sj : a-j ^ Pj.j = \,...,S 

Random grammar includes also positive numbers Aj (rates). That is at time interval (i, t + dt) in the word a{t) 
any subword aj is independently replaced by /3j with probability Xjdt (continuous time Markov chain). 

For graph grammar a{t) are marked graphs, aj,Pj are (small) connected marked graphs. Thus, aj is 
deleted from the graph and 13 j is pasted instead (some restrictions needed of course). Note that ordinary 
grammar is a particular case, corresponding to linear marked graphs. 

One of the problems - invariant measure and conserved characteristics with respect to given graph grammar 
dynamics was considered in \12i |2T1 [22] . 

3 Quantum Graph Grammar 

What is quantum graph [26\. Consider Hilbert space 12{{G}) with (orthonormal) basis ec, enumerated by finite 
graphs G. Or by finite marked graphs if the set of marks is finite. First example is linear marked graphs - 
quantum words. 

To define quantum dynamics assume that if aj — (a^ —>/?;)£ Sub then also inverse substitution a* — (Pi — > 
ai) G Sub. Denote Sj{k) the substitution Sj applied to subword of the word a starting on fc-th symbol of the 
word a. Introduce the Hamiltonian 

\Sub\ oo 

i=i fc=i 

The first simple result is: this Hamiltonian is selfadjoint in 12{{G}), that is the quantum evolution is well-defined. 

Gibbs and Quantum Spaces G - class of finite graphs with a function f : V ^ R , called spin graphs 
{G,f), Gn - class of such spin graphs of radius < D. Potential is defined as some function $ : G^i R. 
Hamiltonian H : G ^ R is 

^((G,/))-E*(^) 
where the sum over all sub spin subgraphs of (G, /). Partition function 

Zn^ J2 exp(-/3ff(G,/)) 
(G,/)eGN 
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Gibbs measure on G/ at 

/iw(G,/) =Z^iexp(-/3iJ(G,/)) 

There are many results-examples (by physicists and mathematicians) related to "quantum gravity", see for 
example [T7t l25j and references therein. 



4 Trading network as Boltzmann mechanics of communicating vessels 

Standard financial mathematics considers games of one or small number of players against the chance (random 
market). Recently, a new approach (called multi-agent models) appeared which considers the games of many 
players against each other. This theory is at the starting point and its models are mainly mean-field models. 

In this section some local models are considered where there are many players and many financial or 
trading instruments. Our model develops simpler models of (|23|, 124)). The model resembles communication 
and transportation networks - the main difference is that the nodes have special dynamical values (moving 
boundaries, or real prices). The clients have also their own subjective prices and their interaction (transaction) 
with the nodes depend on these prices. This model does not describe any real situation (and any other existing 
multi-agent model as well) but we hope that some features of this model will be useful for future more realistic 
models. 



Free one-phase Boltzmann dynamics Consider the phase space S = I x Iq, where / C i? is an infinite 
interval and Iq ~ [— Vb,Vb], < Vq < oo. On S at any time t > a random locally finite configuration 
{{xi{t), Vi{t))} of particles is given with coordinates Xi ^ I and velocities vi G /q. Assume that this configuration 
at any time t has distribution Pt with one-particle correlation function f{x,v,t) defined so that for any subset 
A C S of the phase space 

Ei^{i: {xi,Vi) £ A} ^ / f{x,v,t)dxdv 

J A 

One can have in mind Poisson measure Pq at time t — 0. Any particle moves always with its initial velocity, 
independently of other particles. Also there is Poisson income fiow of particles from exterior with rate X{x, v, t), 
that is during time interval [t, t + dt] the mean number of incoming particles to the cell [x, x -\- dx\ x [v.v + dv\ 
of the phase space is A(2;, t)dxdvdt. Assume moreover that each particle can die (disappear) with exponential 
distribution having rate /i(x, w, t). This means that during time dt ^{x, v, t)dxdvdt particles leave the cell dxdv. 
Remind that we assume boundedness of velocities, that is 

f(x, V, t) = X{x, V, t) = n{x, V, t) = 0, \v\ > Vo 



Lemma 1 For any x G I and t < '^^ y^"^'* ; where d{x,dl) is the distance of the point x from the boundary of I , 
the standard linear Boltzmann equation holds 

% + = ^^(^' ^' ^)/(-^' ^' ^) + ^(^^ *) (1) 
This is trivial for /i = A = 0. In fact, for small 5 > we have 

fix,v,t + S) = fix-v6,v,t) (2) 

if X is not on the boundary of / and 5 is sufficiently small. Subtracting f{x, v, t) from both parts of this equality, 
dividing by 5 and taking the limit (5 — 0, we have 

df df 

The unique solution of the Cauchy problem for ^ is 

f{x,v,t) = f{x - vt,v,0) 

If A 7^ 0, /.t = fi{x, v) ^ then it is also easy to see that the equation ([1]) holds. Note that if A = and fi does 
not depend on t, there is also explicit solution, see section XI. 12 in |27] 

f{x,v,t) — f{x — vt,v,0)exp{— / ^{x — vs,v)ds) 

Jo 
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Two phases - particle dynamics We shall define two types of dynamics - particle dynamics and continuum 
media dynamics. 

In the particle dynamics (±)-phases consist of (±)-particles so that each (~)-particle is to the left of any 
(+)-particle. Denote b{t) £ R (boundary between phases) the coordinate of the leftmost (+)-particle. Then 
for X > h{t) there is (+)-phase and for x < b{t) there is (— )-phase. Particles move, as above, with their own 
velocities until a (— )-minus particle reaches the point b{t), then it disappears together with the (+)-particle at 
b{t) and the point b{t) jumps to the coordinate of the new leftmost (+)-particle. After this, the process proceeds 
similarly. 

Random configurations of particles are defined by the correlation functions f±{x,v,t) correspondingly. As- 
sume that also the functions A-|-(r, i), /i±(r, t), r > 0, are defined, smooth on i?+ and zero if r > Rq for some 
< i?o < oo. 

The dynamics of one point correlation functions f±{x,v,t) for x ^ b{t), that is on (6(i),oo) and (— cxD,6(t)) 
correspondingly, is given by the equations (already non-linear as b{t) is unknown) 

^+1^^ - -^l±{x -b{t),v,t)f±{x,v,t) + \±{x ~b{t),v,t) (4) 

This means that we assume that arrivals and departures depend only on the distance r — \x ~ b{t)\. 

Thus two phases add reactions between particles of different phases. The following interpretation is useful. 
We consider one instrument (stocks, futures, houses or other real estate etc.). There are two types of 
traders - (+)-particles correspond to sellers and (— )-particles to buyers, Xi are subjective prices comfortable 
for the trader i. Collision between particles corresponds to transaction, after this both leave the market. In 
more general cases it will be possible that they do nor leave the market (see below). 

We consider here a particular case when for some constant velocities v± and for any t 

f±{x,v,t) = p±{x,t)S{v - v±) 

For this to hold at any time t it is sufficient to demand that this holds for t — 0. Initial conditions are defined 
by the initial densities p±{r, 0). The velocities v± can be interpreted as averaged velocities for sellers and buyers 
correspondingly. 



Two phases - fluid dynamics It can occur that under some scaling the defined particle dynamics tends to 
some kind of continuous (fluid) picture, see [23], but we shall not pursue this way here. Instead, we consider 
continuous densities of (-f)-masses and (— )-masses and shall define their dynamics directly. We assume that 
at each time t there exists point b{t) - boundary between phases. There are two phases with initial densities 
p+{r,0), p^{r,0) where 

r ^ r{t) = \x- b{t)\ = ±(a; - b{t)) 

correspondingly. Phases move with velocities w-t correspondingly. Collision of plus and minus masses (at the 
point b{t)) leads to their cancellation in equal amount. There is more realistic possibility - to make the cancel- 
lation proportional to the current price, but we do not consider this possibility here. 

We obtain equations for the triple {b{t), p+{r,t), p-{r,t)) similarly to the way how the equations of contin- 
uum mechanuics are derived in the textbooks, that is using conservation laws. Here there is only one - mass 
conservation law. 

First of all, obtain the equation for the boundary. Assume b{t) smooth and put (3 — Then for time dt 

the amount of positive mass, reaching the boundary will be 

M+{l3,t)dt= / f+{r,v,t)dv + o{dt) ^ dt f+{0,v,t){-v + l3)dv + o{dt) 

Jv-p<0 Jr<{-v+P)dt Jv-/3<0 

111 fact, income and outcome give the contribution o{dt). Similarly for negative mass 

M_{(3,t)dt= 1 1 f_{r,v,t)dv + o{dt) =dt 1 f^{0,v,t){v - f3)dv + o{dt) 

Jv~l3>0 Jr<{v-I3)dt Jv-fi>0 

Lemma 2 For any t there exists unique /3 — f3{t) such that 

M+{P,t) ^ M.{p,t) (5) 
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In fact, consider the equation with respect to /? 

f f+iO,v,t){~v + P)dv^ f f-{Q,v,t){v - f3)dv 

Jv-I3<0 Jv-fS>0 

Then if (3 increases, the right-hand side decreases and the left-hand side increases. 
We can rewrite the equation ([S]) in our case 

p+(o, t){~v+ + m) = p- (0, - m) (6) 

from where we can get j3{t) 

P+{0,t) + p^{0,t) 

Now we should write the equations for the densities. For p^{r,t) we get 

p+ir, t + At)^ p+{r - {v+ - m)^t, t) - ^^+{r, t)p+{r, t)At + A+(r, t)At + o{At) = 

= p+{r,t) - {v+ - i3{t))^^^^^At - /i+(r, t)p+{r,t)At + X+(r, t)At + o{At) 
In the limit Ai ^ 

dp+{r,t) dp+{r,t) t n , \ f ^\ fa\ 
Q-^ = I3[t)) g;^ M+ t)P+ {r, t) + \+ [r, t) (8) 

Similarly p_ (r, t) : 

p-{r,t + At) = p_(r + (w_ - (i(t))At,t) ~ p,_{r,t)p+{r,t)At + X_{r,t)At + o{At) 

= iv^-pit))^P^-^^ir,t)p^ir,t) + X^ir,t) (9) 

It would be nice to prove accurately that the solution of equations (|6l8l9p exists for any t > and is unique, 
but we did not try to do this. 

Fixed points and stationary points Assume that the functions A±(r) = X±{r,t) and n±{r) — p.±{r,t) do 
not depend on t (remind that they were assumed to have compact support). Denote 

7ir!"' = -«+^^ A+(x)exp^-^y /i+(y)dy^ dx, 7^7' = wZ^ y A_(x) exp ^ p.^{y)dy^ dx 

and 

( (+) ^'-7crM 

7cr = max -I'^J' , 

V / 

We define the fixed point of our dynamics by the conditions: /3{t) = and p±{r,t) do not depend on time. 
Alternatively the fixed points are defined as any solutions of the stationary version 

p^.(0)^;++p_(0)^;_ -0 (10) 
- ^+^^ - M+Mp+(r) + A+(r) = (11) 

V-^^^-p-{r)p-{r) + \^{r) = Q (12) 

of the system (|6l8l9p . We will prove that there exists a family of fixed points depending on a real parameter. 

Similarly, we call stationary point any solution of the system of equations (|6I8I9I) . where /3 — (3{t) and the 
densities do not depend on t. We shall prove that there is a family of stationary points depending on two real 
parameters. 

We say that a fixed (or stationary) point has finite mass if 

/•OO 

/ p±{r)dr < 00 
Jo 
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Theorem 3 Let the parameters X±{r), fj,±(r) and v±be fixed. Then 

1. For any value of the parameter 7+ — P+(0) there is at most one fixed point. For 7+ < there is no any 
fixed point. For 7-|_ > ^cr there exists exactly one fixed point defined by 

p+{r) ^ exp(^--^J fi+{x)dx^ (^p+{0) + v^^ X+{x) exp (^-^ j p+{y)dy^ dx^ (13) 

P-{r) = wZ^exp^-!-^ p^(x)dx^ (^-v+p+{0) - j X^{x) exp (^--^ J fi^{y)dy^ dx^ (14) 

2. The fixed point has finite mass if jct^ — jcr ^ 

3. For any ^^,^^such that 

7+ = P+(0) > 7i+\ 7- = P-(0) > 



there is exactly one stationary point. Then the densities are defined by formulas j^[ ) and the boundary 
velocity is 

P+(0)«++P-(0)«- 
p+(0)+p-(0) 

4. Stationary point has finite mass iff J+ = Jcr \ 7- — Icr ^ ■ 

Proof. Solving equations (|8l9p we get for any r > equations and Note that, by equations 

and ([T^ . densities p-{r), p+{r) are positive iff 7+ > 7- > ^if'' ■ Taking into account equation ([TU]) we 

get the first asserion of the theorem. 

For the stationary points the densities are again defined by equations and (|14p . We have two conditions 
for them to be non- negative. Then the boundary will move with constant velocity defined from equation ((T]). 

More complicated one market model Note that collision of masses of two phases create total annihilation 
flow 

= («_ - m)P-{^.i) = - /3(i))p+(0,i) 

of the disappearing (±)-particles. Here we assume that a part of annihilating particles does not disappear 
but can transform to particles of the other phase jumping from the collision point to some point r. On the 
language of continuous media this means that there are output flows of mass v{+, — , r, t) and +, r, t) such 
that 



p(+,-,r,Odr < 1, / p{-,+,r,t)dr <l 
Jo Jo 

where 

-,r,t)^ ^^(ii-fl^, Pi-, +, t) . ^-t^ 
For such model we have the system of three equations 
t;_p_(0,0+w+P+(0,0 



m = 



p-{0,t) + p+{0,t) 



^^^^ = - {v+ -m) ^^^qI' M+ {r, t)p+ (r, t) + A+ (r, t) + (^;_ - /3(t))p_ (0, t)p{-, +,r,t) (15) 

^^^^ = (v^ - p{t))^P^^ - p^ [r, t)p^ (r, t) + A_ (r, t) - {v+ - l3{t))P+ (0, t)p{+, -,r,t) 

We again assume that the functions p±{r,t), X±{r,t),p{—, +,r,t),p{+, —,r,t) do not depend on t and have 
compact support. Introduce the functions 

1 r 1 r 

F+{x) = / p+{y)dy, F^{x) = — / p-{y)dy 

v+ Jo v_ Jf) 

Denote 
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p{—,+,x)exj>{—F^{x))dx, Q!-| = / p(+, —, x) exp {—F^(x)) dx 

Jo 

and assume that a a-| < 1. Define 

/ (+) (-) \ 

7cr = max 



1 — a — [- ' — — a-i — ) / 

Theorem 4 Let the parameters X±{r), p±{r), p{~,+,r,t),p{+, —,r,t) and v±be given. Then 

1. For any value of the parameter 7+ = /9+(0) there is at most one fixed point. For 7+ < ^^r there is no any 
fixed point. For 7+ > ^cr there exists exactly one fixed point. It is 

p+{r)=c^p{F+{r)) (^p+{0) + v^^ {\+{x) - v+p+{0)p{-, +, x)) cxp {-F+{x)) dx 



P-{r) = vZ' exp (F_(r)) ^-v+p+{0) ~ (A_(a;) - v+p+{0)p{+, -,x)) exp {-F^{x)) dx 
2. There is a unique stationary point with finite mass. It is 

p+{r)^exp{F+{r)) (^p+{0) + {X+{x) + {v^ - I3)p^{0)p{-, +, x)) exp {-F+{x)) dx 

p^{r)^exp{F^{r)) (^p^{0) - vZ' {X-{x) - {v+ ^ p)p+{0)p{+, -,x)) exp {-F.{x)) dx 
and we denote 



— — a — ^) — fia — I 
(-) 



P-(0) 



where (3 is a root (belonging to the interval of quadratic equation I12S\) . It exists and is unique. 

Proof. 1. Similarly to the first part of theorem [31 

2. As follows from system (|15p the equations for the stationary points are 

= {v--P)p-iO) + {v+-P)p+{0) 

^ _„^^^_^^(^)p^(^) + A+(r) + («_-/?)p_(OM-,+,r) (16) 
or 

= v^^P-^-^_{r)p^{r)+X^{r)-{v+^(3)p+{0)p{+,-,r) 
Solving these linear first order equations we get 

p+(r)=exp(F+(r)) (^p+{0) + v^' {X+{x) + {v_ - l3)p.{0)p{-,+,x)) exp {-F+{x)) dx^ (17) 

p_(r) =exp(F_(r)) (^p^{0) - vZ' (A_(x) - {v+ - f3)p+{0)p{+, x)) exp {-F^{x)) dx^ (18) 
We are looking for a stationary point with finite mass such that 

p±{r)dr < 00 (19) 
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Then by (|T7| . ([T8| . (fT9| . |T6]), a stationary point is uniquely defined by three parameters 7± — p±(0), /3 which 
satisfy the following equations 



= ^-7^7' - («+ - /^)7+a+- (20) 
(w- - /3)7- = - /3)7+ 

where /? S (w+ , W- ) . We show that this system has a unique solution. Using the third equation of the system, 
we get from the first two 

(+) 

^+ " T^ A « V^^) 

7 = (22) 

Substituting these expressions to the third one we come to the quadratic equation with respect to 

(CT+a+_ - (T_a_+)(— u+ + fi){v^ - ^) + — cr+V-)^ — u+u_(cr_ - (7+) = (23) 

where, for shortness, we denote cr+ = — u+7cr^\ (7_ = v-jir ^■ 

Consider first the case when cr+a-f — a-a [. . Note that the boundary velocity should satisfy 

v+ < l3 < v^. One can show easily that there is always one root of the equation in the interval < (3 < v-. 
Now one should verify that 7+ , 7_ , defined by (PT|) and are non-negative. By (PT|) (1^^ one of the values 
7+,7_ is always positive. Then by the third equation of the system (I20p also the other value is positive as 
V- - (3, -v+ + /3 > 0. Thus there exists the unique fixed point satisfying ((17]), ([18]), ([21]) (|22|) . 

Is cr+Q!-i — (7_a y = 0, we have a linear equation with respect to /?, we gives 

/3 = 



and from (|2ip ([221) we get 7+ = cr+f.,.^ and 7_ = a-v_^. In this case also a stationary point exists and is 
unique. 

Networks with many markets Let us call the previous model an elementary market. A network is a set V 
of elementary markets with similar parameters and variables indexed hy m ^ V 

{r,t),b^{t),f3ra{t) 

There are also other parameters interconnecting the markets. Denote z^+.„i(i) {v-,m{t)) the total annihilation 
flow of (±)-particles from the market m. As they are equal we denote Vmit) = ^+,'m{t) = ^-,m{t)- Let 

Vk,mi + , +, r, t), Vk,rn{ + , r, t),Vk,m{~, +, T, t) , Vk,m{~, +, r, t) 

be the parts of these annihilation flows of (±)-particles, that after the transaction on the market m, become 
(=F) -particles on the market k with the coordinate r. Denote 

Pkm(±,±,r,t) = 

We mean that pkm{+, +, t) — Pkmi~j ~i ''i = 0- Then for any k and t the conditions 



J2 / Pkr,r{+,-,r,t)dr<l, / 



Pkm{-,+,r,t)dr < 1 
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should hold. Denote by \V\ the cardinality of the set V. We have then the following system of 3|y| equations: 

^— ,mP— ,m -\- ,m P -\- .m 

{0,t) 



dp+^ra{r,t) dp+^m{r,t) f , \ ( ^\ 

= -(w+,m - Pm(i)) g;^ fi+,rn{r,t)p+^rn[r,t) + \+^m,[r,t) 

+ -/3fe(t))p_,fc(0,tK™(-,+,r,t) 

dp-^„,{r,t) dp_^^{r,t) / +\ , \ r 

ot or 

- ^{v+,k - l3k{t))p+,k{Q,t)pk,n{+,-,r,t) 

k£V 

Fixed points Again we assume X±^mir, t), iJ,±^m{,r, t),pkmii, ±, r, t) do not depend on t and have a com- 
pact support. Put 

Jo Jo 
X+,m^J X+^rn{x)e^p(^-F^\x)jdx, X^^rn^ J A_ (x) cxp (^-fI^^x) j (24) 

afcm(-,+)= / Pkrn{-,+,x)c^p [-F!f'\x)j dx, akni{ + ,-) = / Pfcm ( + , - , x) CXp ^-FI™-* (x) j 



, , , , , , . - , , dx 

Jo ' ' \ ^ " ' J 'Jo ' ^ ~ 

for fc, m G 

Define matrices A A_| with elements ckfcm(+j~)! where k,m ^ V . and assume, that they 

have the following property: 

Vfc "fc™(±'±) < 1. 3fco ^ afc™(±,±) < 1 (25) 

For two vectors a = (a^) and 6 = (6^) we shall write a > b{a > b) if a.; > hi [a^ > 6^) for all coordinates. Consider 
the following system of inequalities with respect s 

s{E - A^+)>X+,s{E - A+^)>X- (26) 

where E is the identity matrix and A± are vectors with coordinates X±^m defined by (P^ . We say that this 
system has a positive solution if there is vector s with positive coordinates satisfying both inequalities in ([26| . 

Generally, this system may not have a positive solution. If one of the matrices A |_, A^ is diagonal or zero 

the set of positive solutions is nonempty. 

Theorem 5 Each solution s ~ (sm, m Cz V) > of the system ^26\} uniquely defines the fixed point as follows: 
p+,m{r) = -u+^^„exp (^F|'"-'(r)^ ^s™ - J ^A+,™(a;) + ^ SkPkmi-,+,x)^ exp (^-F'i^\x)j dxj 

p^,m{r) = vZ^m (-^'-'"V)) - J ^A_,,„(x) + ^ SkPkm{+, ~, x)^ cxp (^-F["^\x)j dxj 

If the set of positive solutions of system I126\) is empty there is no any fixed point. 
Proof. The fixed points satisfy the system consisting of 3|1^| equation: 

= W-,mP-,m(0,t) + (0,0 

= -V+^m ^^^^^^^ - fl+,rnir)p+,mir) + X+^rnjr) + t'-.fcP- ,fc (O)pfcm (~ , + , r) (27) 

kev 

dp (t*) 

= V-^jn ^ fJ.-^rn{r)p-,m{r) + A-,m(r) - V+^kP+,k{^)Pkm{+, -, t) 

kev 
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Solving first order linear differential equations we get 



P+,mW = exp (F|"'(r)) I p+,r„(0) + v^^^^ / | A+^,„(x) + ^ t;_,fcp_,fc(0)pfc„(-, +, x) j exp (^-F^"'\x)j 



kev 



dx j 
(28) 

P-,m(r-) = exp (Fl™'(r)) | p_,m(0) - vZl,n / f A_,,„(a;) - ^ w+,fcp+,/c(0)pfc™(+, j exp (^-F^"'\x)j 



kev 



dx j 
(29) 

for m V. 

Using equations = w-.-m/O-.m(0) + w+,mP+,m(0), we conclude that solutions (l28l) . (l29l) are uniquely defined 
by parameters = — u_|_.j„/3+^m(0), m G V^, and one can write 



p+,m(r) = -w+^^^exp (^F|'"'(r)^ y^m^J |^A+,m(x) + SfcPfcm(-, +, a^) j exp (-F^'\x)) dxj (30) 
p_,„(r) = vZl,r^ exp (fI^V)) (^s™ - ^ |^A_,™(a:) + SfcPfc„,(+, a:)^ exp (-f1™^(x)) d.T^ (31) 



Whereas the densities (|30p . ([31]) are nonnegative for all r > the following conditions must be satisfied 

(-,+) (32) 



Sm > A-,TO + Skakm{ + , -) (33) 



for all m G V. These inequalities are equvalent to system ([26]). 

So the fixed points exist iff there exist positive solutions of system ([26| . 
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